


AY 1) WFEOROLOGICAL 
MPN AZINE 


HER MAJESTY’S April 1979 
STATIONERY 
OFFICE Met.O. 922 No. 1281 Vol. 108 








METEOROLOGICAL OFFICE Met.0.922 


THE METEOROLOGICAL 
MAGAZINE 


No. 1281, April 1979, Vol. 108 





$51.577.33: 519.24].25 


A spectral and filter analysis of long-period rainfall records in England and Wales 


By R. C. Tabony 
(Meteorological Office, Bracknell) 


Summary 


Four long-period rainfall records in England and Wales have been subjected to maximum entropy spectral analysis 
in order to identify any regular periodicities which may be present. Unitary filters were used as a coarse form of spectral 
analysis in support of the maximum entropy technique, and the series were filtered to examine the temporal variations 
in the amplitudes of some of the periodicities. A wide range of seasons and epochs was examined, and although some 
regular periodicities were evident, random fluctuations were dominant. Cycles that could be regarded as permanent 
were located at 2-1 and 2-4 years, of which the latter, associated with the summer half-year, was the most important. 
Cycles that were less prominent, but which attained 5 per cent significance over the complete length of all the series 
examined, were found at periods of 3-9, 5, and 6 years. 


1. Introduction 


Major searches for periodicities in meteororological data were carried out in the 1920s and 1930s, a 
notable investigation being due to Brunt (1927). Later the search for cyclic behaviour fell from favour, 
but the advent of computer techniques renewed interest, and Ward and Shapiro (1961) made power 
spectral analyses of a wide variety of meteorological parameters on a world-wide basis. A general 
review of investigations into cyclic behaviour in the atmosphere is provided by Lamb (1972), but the 
majority of studies have dealt with temperature or pressure, Recently, the well-known rainfall record 
for England and Wales has been augmented by three homogeneous series representative of particular 
sites in England. In addition, the last 20 years or so have seen rapid developments in the field of power 
spectral analysis, with the latest techniques providing much finer resolution than the earlier methods. In 
the search for periodicities in rainfall in this paper, therefore, the latest techniques of spectral analysis 
are combined with the long-period rainfall series which have only recently been published. 

In this paper the complete lengths of the four rainfall series are spectrally analysed and the significant 
spectral peaks are tabulated for a number of seasons. The analysis was repeated on component 80 year 
epochs of the series to determine the degree of permanence of the spectral peaks. A quasi-orthogonal 
unitary filter analysis applied to the complete records was used as a coarse form of spectral analysis in 
support of the maximum entropy technique, while filtered time series were used to examine in more 
detail the amplitude of some of the periodicities found. 
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2. Data 


The four series of monthly rainfall totals analysed were 

(i) England and Wales from 1727, first compiled by Nicholas and Glasspoole (1932) and since 

maintained by the Meteorological Office. 

(ii) Kew from 1697, first prepared by Wales-Smith (1971) and revised in 1978. 

(iii) Pode Hole, Lincolnshire, from 1726, assembled by Craddock and Wales-Smith (1977). 

(iv) Manchester from 1786, published by Manley (1973) and revised in 1976. 

Decadal means of annual rainfall for the four series are presented in Figure 1. The two main features 
of interest are the rise in the totals for England and Wales, most of which occurred before 1820, and the 
convex nature of the graph for Pode Hole, with low values in the early eighteenth century and high totals 
in the first half of the nineteenth century. These features contrast with the relative absence of long- 
period trends at Manchester and Kew. However, the relatively large changes in rainfall indicated by the 
series for Pode Hole and England and Wales may not be real. The first 40 years of the England and 
Wales series are based mostly on data from only two or three stations with unorthodox sites by modern 
standards. During the next 60 years the number of stations used was only about 10, and they are not 
evenly distributed over the country. From 1820 onwards, the number of stations was sufficient to enable 
a mapping technique to be used to derive the areal rainfall, and thus it is only from this date that the 
England and Wales series can be regarded as reliable. The compilation of the series for Pode Hole, a 
rural site, was inevitably based on less evidence than the series for Kew and Manchester, yet both the 
latter have been subject to revision. It may well be that the Pode Hole series will require revision in 
future, with the probability that the large fluctuation in decadal rainfall in the present record will be 
diminished. 

The uncertainty attached to the values accorded to the earlier portions of the records casts doubt on 
the reality of trends or long-period cycles in the data. Fluctuations of shorter periods will be less 
affected, however, and cycles of wavelengths less than about 25 years may be regarded as being 
relatively undistorted. 


3. Spectral techniques used 


3.1 Brief survey of methods of spectral analysis. When a series of observations is transformed from 
the time to the frequency domain, and the square of the amplitude (i.e. power) is plotted against fre- 
quency, the result is known as a power spectrum. The area under the curve in a power spectrum is pro- 
portional to the variance. A good account of spectral analysis in general is provided by Bath(1974). The 
transformation from the time to the frequency domain may be made by taking the Fourier transform of 
the time series (direct method) or of the autocorrelation function of the series (indirect method). 

In the pre-computer era, harmonic analysis (the direct approach) was used, and the results were 
usually displayed as a plot of amplitude against frequency (the periodogram). Later, Blackman and 
Tukey (1959) developed a method for computers based on the indirect approach, but Cooley and Tukey 
(1965) then derived the Fast Fourier Transform (FFT). This is essentially the direct method, but the 
ordinary ‘direct’ formulae are replaced by more computationally efficient ones. These methods are 
efficient when operating only on M* data points where M can equal 2, 3, 4,.... In practice, therefore, 
the series is often truncated or extended at both ends with mean values until the number of points equals 
a multiple of M. 

Jenkinson (1977) has described a direct method which produces a quasi-continuous form of the 
periodogram. The increase in resolution obtained is especially pronounced at long wavelengths because 
the frequency elements are chosen to be logarithmically spaced. 
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Figure 1. Decadal means of annual rainfall. 
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The maximum entropy method, first suggested by Burg (1967), is described by Lacoss (1971) and 
Ulrych and Bishop (1975). It is an indirect method, based on the transformation of an autoregressive 
process which neither adds information to nor subtracts information from the data (the concept of 
maximum entropy). The resulting extension of the original data series means that the method is capable 
of much higher resolution than the earlier methods. 

In this paper a maximum entropy package due to Ross (1975) has been used. 

3.2 Calculation of variance and significance of spectral peaks. Spectral peaks have arbitrarily been taken 
to lie between those frequencies where the power is one-third of its peak value. The variance attributed 
to a given peak has therefore been limited to the area contained between the one-third power points, as 
illustrated in Figure 2. 


Y Variance between 
GY one-third power 
points 
Rex] Variance between 
22S One-third power 
points expected from 
a random series 








—> Frequency 





Figure 2. Definition of spectral peak by one-third power points. 


The estimation of the statistical significance of a peak is based on the variance ratio, that is to say the 
ratio of the variance observed to that expected from a random series. If the width of the peak defined by 
the one-third power points is Af, the mean power over the frequency range Afis h,, and the mean power 
over the whole of the spectrum is h,, then the variance ratio R is given byh,/h,. As serial correlation 
was absent in the rainfall data examined, the number of degrees of freedom associated with the frequency 
range 0 to 0.5 was set at N, the number of data points in the original series. The degrees of freedom were 
assumed to be uniformly distributed in frequency space, so the number attributed to the frequency range 
Af was taken as n = N x Af/0.5. The significance of a peak was estimated from statistical tables of 
Snedecor’s F distribution using the variance ratio R and the number of degrees of freedom N and n. 

There are two important criticisms of the procedure described above. The first is that the use of the 
F test to assess the significance of the variance ratio R applies only when the frequency range Af is 
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arbitrarily defined as, for example, when the frequency range 0 to 0.5 is divided into M equal parts. By 
choosing the width and central frequency of Af to suit the peak under examination, a systematic over- 
estimation of significance results. In this paper, however, the significance of a given peak relative to 
another is just as important as its absolute value, and the procedure described above leads to more 
accurate values of this than if the spectral peaks had been assessed over arbitrarily pre-defined frequency 
ranges. 

The second criticism relates to the uncertainty in the number of degrees of freedom when N x Af/0.5 
is close to unity. This uncertainty only applies to narrow peaks for which Af is small, however, and in 
most cases dealt with in this paper the estimate of N x Af/0.5 will be reasonable. 

3.3 Comparison between maximum entropy and other methods of spectral analysis. The maximum 
entropy method of spectral analysis is not universally accepted, and so itis pertinent to make comparisons 
with other methods, such as the FFT and Jenkinson periodogram. In the maximum entropy method the 
number of autoregressive coefficients used to extend the data has to be chosen by the user, and this is 
perhaps the most frequent criticism of the technique. The number of coefficients determines the resolu- 
tion or ‘spikiness’ of the spectrum. A small number of coefficients produces a highly smoothed estimate 
but the resolution becomes finer as the number of coefficients is increased, as illustrated in Figure 3. 
With other methods the resolution of a spectrum increases with the number of data points N, and is fixed 
for a given length of input data. Akaike (1970) suggested that the number of coefficients to be used in 
the maximum entropy method should be such that any further increase would produce a reduction in 
variance that is not statistically significant. This is of little value, as it determines the length of the auto- 
regressive series from a measure of its predictive power. A much greater length may be needed to dis- 
play periodicities of scientific interest which may not have much predictive power. A practical solution 
is to use N/3 coefficients for N less than 100, and then to use a decreasing fraction of N as N increases 
beyond 100. 

The significant spectral peaks obtained from an analysis of annual rainfall at Kew from 1697 to 1975 
are presented in Table I. The results are obtained from the FFT, Jenkinson periodogram, and the 
maximum entropy method using 25, 40, 60 and 100 coefficients. The variances and significances were 
obtained using the procedures described in section 3.2, and only peaks which were significant at 5 per 
cent or better have been entered in the table. 

Considering first the results of the maximum entropy method, the number of peaks produced in each 
spectrum was about 40 per cent of the number of coefficients used (see Figure 3). Thus, as the number of 
coefficients and the number of peaks increases, the proportion of the variance accounted for by each 
peak decreases. In contrast, the estimated significance of a given peak increases with the number of 
coefficients. This fact, combined with the greater number of peaks produced for the larger number of 
coefficients, means that the number of significant peaks produced increases rapidly with the number of 
coefficients. Thus a tabulation of the significant spectral peaks derived from a maximum entropy 
spectral analysis is highly dependent on the number of coefficients used. A large number of coefficients 
produces a relatively large number of significant peaks each of which accounts for only a small propor- 
tion of the variance, while a small number of coefficients produces a relatively small number of signifi- 
cant peaks each of which accounts for a more substantial portion of the variance. Despite this variation 
in the absolute values of the variances and significances associated with a given fluctuation, the relative 
importance of the peaks is independent of the number of coefficients. In the analysis of long-period 
records which follows, the maximum entropy method has been used with 26 coefficients for 80 year 
epochs, 60 coefficients over the complete length of record (about 250 points) and 100 coefficients for the 
complete records divided into three month seasons (about 1000 points). In all cases, the linear trend 
was removed from the data before the power spectra were produced. 
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Figure 3. Maximum entropy power spectra of annual (calendar year) rainfall at Kew from 1697 to 1975 using 100, 60, 
40, and 25 coefficients. 
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Table I. Significant spectral peaks of annual rainfall at Kew from 1697 to 1975. Comparison of FFT, 
Jenkinson and maximum entropy methods of spectral analysis 


FAST FOURIER eer 
> 6°02 3-91 93 
3 4:8 5:2 3. 1 
5 5 3 


JENKINSON PERIODOGRAM 
m 3-90 2°95 2: 
= 4-7 2°3 2 
5 + a 
MAXIMUM ENTROPY (100 COEFFICIENTS) 
12-3 6:00 3-89 2-94 2°85 46 
5 


z 
3-8 2-6 4:5 2-2 2:2 2: 

3 a a = 5 3 
MAXIMUM aa (60 COEFFICIENTS) 
12-6 6°11 

78 3°4 rf I 

1 5 5 
MAXIMUM ENTROPY (40 COEFFICIENTS) 
12-3 


6:0 
5 


MAXIMUM ENTROPY (25 COEFFICIENTS) 


1 


1 
9- 
5 


46 
- 


T = period (years); V = proportion of variance (%); S = level of significance (%). Peaks are listed only if significant 
at 5 % level or better (similar remarks apply to Tables V- XIII). The horizontal displacement from the left-hand margin 
of a column in this and similar tables is roughly proportional to the associated frequency, i.e. inversely proportional to 
the period. 


The results of the FFT and the Jenkinson periodogram agree well with one another and also with the 
maximum entropy analysis for 60 coefficients. Table I shows that an analysis made by any of these 
methods will enable the variance associated with a given fluctuation to be estimated to within a factor of 
about 2, and its statistical significance assessed to within a factor of about 5. The degree of absolute 
accuracy is therefore poor, but the assessment of the relative importance of peaks is much better. The 
main aim of this paper is to analyse rainfall data to see if any periodicities of a permanent nature are 
present and if so, whether they account for a useful proportion of the variance. The results presented in 
Table I show that if any major periodicities are present, then all three methods of spectral analysis 
examined here will detect them, and they will also agree as to their relative importance. 


4. Quasi-orthogonal filter analysis 


Sixth-order unitary filters due to Craddock (1968) are used in this section, and the combinations 
chosen are described in Table II (using the notation of Craddock). In a random series, the band pass 
filters labelled 1 to 5 pass approximately 1/6th of the total variance around the peak periods indicated in 
Table IL, while the low and high pass filters labelled 0 and 6 pass approximately 1/12th of the variance 
around wavelengths of oo and 2 respectively. The variance obtained from a given series can be compared 
with the variance expected to be obtained from a random series, and the ratio between the two is known 
as the normalized variance. The significance of a departure of the ratio from unity is assessed using 
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Table I. Description of quasi-orthogonal filters used 


Number of Peak period Half-power points 
filter Composition Number of data points 
0 1:0 Fy, ; . 23-5 
S 18-6 7:50 
7-88 4-74 
4-81 3-42 
3-49 2-68 
2-72 2°24 
2°19 


209999 
NNNNNN 
Pa Pe Be | 
2eana 2 @ 


Snedecor’s F test, as first suggested by Craddock (1957). The results of an analysis of rainfall for 
England and Wales, Kew, Pode Hole, and Manchester for a variety of seasons are presented in Table III. 
In this paper, the word ‘season’ is used to describe a period of consecutive months starting in a specified 
month. The summer half-year refers to the period April-September and the winter half-year to October— 
March. The rainfall totals were not serially correlated, and so the number of degrees of freedom 
allocated to each series was equal to the number of data points involved. Only significances of 5 per cent 
or better have been entered in the table. 

In Table II most of the normalized variances are close to unity, showing that the variances of the 
filtered series are not significantly different from those which would be expected from a random series. 
The exception is for filter 0, which measures the variance associated with long-term fluctuations. It can 
be seen that large departures from a random series have occurred at Pode Hole, especially in the autumn, 
and also over England and Wales during the winter half-year. These findings quantify the magnitude of 
the long-period variations in the data illustrated in Figure 1 and discussed in section 2. Apart from those 
with long periods, the only other filtered series to attain significance at the 5 per cent level is that for 
periods of around 2.5 years in the spring at Kew. No great importance is attached to this, as 1 in 20 
random filtered series could be expected to reach the 5 per cent significance level by chance. 

It is interesting to compare the results for 12 months starting in January with 12 months starting in 
July. For longer periods, the behaviour of rainfall totals summed over 12 months will be independent 
of the starting month, but for shorter wavelengths this is not necessarily the case. The biggest differences 
between starting months are found at Kew, where the variance for filters 5 and 6 is less for 12 months 
starting in July than for 12 months starting in January. 

So far only the rainfall occurring in individual seasons has been examined. A chronologically 
ordered data set comprising all seasons mixed together, i.e. spring followed by summer, autumn, winter, 
spring, and so on, would not necessarily give the same results. Accordingly, data sets were formed of 
rainfall summed over the conventional three month seasons with departures from average being ex- 
pressed in terms of the standard deviation for each season. The number of data points was four times as 
many as for the previous series examined, and the shortest period that could be examined fell from 2 
years to 0.5 year. The results of the filter analysis are presented in Table IV, but they show only a 
continuation of normalized variances close to unity. As a result of the half-power points of filter 0 
dropping from 23.5 years to 5.9 years, the normalized variances associated with filter 0 are less than those 
displayed in Table III, and only the series for Pode Hole achieves significance at better than 5 per cent. 


5. Maximum entropy spectral analysis 


5.1 Analysis of individual seasons. A maximum entropy spectral analysis of the complete series for 
England and Wales, Kew, Pode Hole, and Manchester was made for a wide range of seasons—the 
winter and summer half-years and the conventional three month seasons are supplemented by 12 
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Table Ill. Normalized* variance (Vx) and corresponding level of significance (S per cent)} associated with 
quasi-orthogonal filter analysis of seasonal rainfall 


England and Wales 1727-1975 
Number Peak Degrees 12 months 12 months Winter Summer 
of period of from Jan from July half-year half-year Winter Spring Summer Autumn 
filter (years) freedom Vx Vx Vs S Ss Vs S Vy S 
22-0 , ° 2:12 O1 1- 1-80 1 . 1-09 
10- 42-1 P " 1-04 1-12 1-07 
. 43-4 ° . 0-88 ° 0-86 . 1-12 
43-4 , ° 0-95 } 0-80 0-79 
43-4 . 0-83 , 0-86 , 0-88 
42-1 , : 0-76 . 0-91 . 1-30 
22-0 , , 1-08 , 1-11 , 0-87 


Kew 1697-1975 
Peak Degrees 12 months 12 months Winter Summer 
period of from Jan from July half-year half-year Winter Spring Summer 
(years) freedom Vy S Wn +S Vz S§ vy Ss Me S Va S Wa S 
248 0-81 0-88 0-79 0-98 . 0-86 0-72 
47-4 1-04 1-06 1-26 1-06 D 1-01 1-10 
48-8 0-92 1-1 90 1-22 } 1-23 1-03 
48-8 0-89 1-0 0-87 + 0-68 0-82 
1-0 
8 
7 


1 e 

03 
48-8 1-13 “88 0-86 ‘ 0-87 0-84 
47-4 = 1-13 0- “89 1-02 1-35 5 1:24 
248 0-97 0- ‘29 0-97 ‘ 0-92 1:19 


Pode Hole 1726-1975 

Peak Degrees 12 months 12 months Winter Summer 
period of from Jan from July half-year half-year Spring Summer Autumn 
(years) freedom Vy S Ve SS Wa S Vy S Vs S 
22-2 . ‘2 2:01 O2 1-85 1 1-66 4 2:16 0-1 

10- 42:5 : 0-8 0-85 1-19 74 

. 43-7 : , 0-86 1-03 
43-7 } } 1-10 0-96 
43-7 . ° 0-99 0-80 
0-99 


PYYNSB 
oucow~ 


3 0 
1 1 
7 0 
7 0 
7 1 


42-5 ‘ x 0-89 


0- 
0- 
0: 
0- 
22:2 0-83 0-79 0-76 0- 


1- 
1- 
1- 
1- 
0- 
1- 
0- 


Manchester 1786-1975 
Peak Degrees 12 months 12 months Winter Summer 

period of from Jan from July half-year half-year Spring 
(years) freedom Vy S Vx S Ve =S Vai S Wa Siva 3 

16-9 1:65 5 1-72 5 1-11 1-59 0-84 1-57 

10- 323 0-84 0-75 0-94 0-92 1-05 0-96 

: 33-2 0-81 1-05 0-98 0-98 1-11 0-76 

33-2 0-95 1-16 0-82 0-82 0-94 1-06 

33-2 0-93 1-03 1-14 0-69 0-94 0-76 

32:3 0-91 0-68 1-04 1-14 0-94 1-34 

169 0-98 0-72 1-17 1-00 1-09 0-69 


* See section 4 of text. ¢ Listed only if 5% or better. 


months starting in January, April, July and October. The number of coefficients used was 60 and the 
significant spectral peaks are presented in Tables V-VIII. Only peaks significant at 5 per cent or better 
have been entered in the tables. 

Of the features common to all four series, the most impressive are the quasi-biennial peaks at periods 
of around 2.1 and 2.4 years. The most important of these is the 2.4 year cycle which is essentially a 
phenomenon of the summer half-year. This is illustrated in Figure 4, where the large quasi-biennial 
peaks for rainfall at Kew summed over 12 months starting in January (which preserves the summer 
half-year) are seen to disappear when the rainfall is summed over 12 months starting in July (which 
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Table IV. Normalized* variance (Vj) and corresponding level of significance (S per cent){ associated with 
quasi-orthogonal filter analysis of rainfall summed over consecutive three-month periods 


England & Wales Kew Pode Hole Manchester 
1727-1975 1697-1975 1726-1975 1786-1975 
Number Peak Degrees Degrees _ Degrees Degrees 
o period of of of of 
filter (years) freedom Vx S freedom Vy S freedom Vy S freedom Vy S 
©O 88 1-06 99 ; 89 ° 0-6 67 
, 169 1-03 190 . 170 ° 129 
174 0-83 195 175 i 133 
174 1-17 195 | 175 133 
174 0-86 195 i 175 , 133 
169 1-05 190 : 170 0-87 129 
88 0-96 99 , 89 0-92 67 


* See section 4 of text. ¢ Listed only if 5% or better. 
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Figure 4. Maximum entropy power spectra (60 coefficients) of rainfall at Kew summed over (a) 12 months starting in 
January, and (b) 12 months starting in July, during the epoch 1697-1975. 
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Table V. Significant spectral peaks of England and Wales rainfall from 1727 to 1975 
Results were obtained from the Maximum Entropy Method using 60 coefficients. 
T = period (years); V = proportion of variance (%); S = level of significance (%). 
Trend 
12 months starting in January 
3-90 


5: , 4:2 
5 


ezy 
Wwaw 
\o 


12 months starting in April 
3-89 


4:3 
3 


12 months starting in July 
‘92-387 
6 43 
0-9 
12 months starting in October 
95 3-38 
-0 3-8 
0-4 
Winter half-year 
97 
7 


ety 
cot & 
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PsA 
Nv 


taney 


Summer half-year 
3-78 
5-1 
2 
Winter 
3-88 
2:8 
5 
Spring 


Summer 
3-7 


4-7 
5 


divides the summer half-year). The 2.4 year periodicity has previously been identified by Alter (1927) in 
northern Europe and the Pacific coast of the U.S.A., and by Jenkinson (1975) in East Africa. The rela- 
tionship of the 2.1 and 2.4 year cycles to the quasi-biennial oscillation in tropical stratospheric winds is 
uncertain, as there are insufficient data available for the latter phenomenon to enable proper comparisons 
to be made. Since observations became available in 1954, the tropical stratospheric winds have dis- 
played a mean periodicity of around 2.2 years. 

Other periodicities of note include a 3.9 year cycle in annual and winter rainfall, most evident in the 
record for England and Wales, but also visible in the other series. Rainfall summed over the winter 
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Table VI. Significant spectral peaks of Kew rainfall from 1697 to 1975 


Results were obtained from the Maximum Entropy Method using 60 coefficients. 
T = period (years); V = proportion of variance (%); S = level of significance ( %). 


Trend 


12 months starting in January 
3-90 


. oF | 1 
0-1 ° “4 4:1 
as 5 
12 months starting in April 
3-92 : 


Rp 8-6 
5 5 


12 months starting in July 
3-36 


5-0 
4 


12 months starting in October 
3-37 
3-6 
2 


Winter half-year 


“92 
“0 


Summer half-year 
Winter 
Spring 


Summer 


half-year displays a periodicity of 5 years in all the series, but is least well developed at Kew. A period- 
icity around 6 years in annual and summer half-year rainfall is evident at Kew, and although absent 
at Pode Hole is also present for England and Wales and Manchester. 

Peaks important only in individual series include a 50 year cycle in annual rainfall over England and 
Wales (using data from 1727), but this feature is probably detected as a periodicity of 40 years at 
Manchester (using data from 1786). A 2.9 year cycle is also prominent at Manchester for annual and 
winter rainfall, while the Kew data indicate a 12 year periodicity in annual and summer rainfall. At 
Pode Hole peaks at 3.3 years and very long periods are present, the latter being the attempt by the 
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Tabie VII. Significant spectral peaks of Pode Hole rainfall from 1726 to 1975 


Results were obtained from the Maximum Entropy Method using 60 coefficients. 
T = period (years); V = proportion of variance (%); S = level of significance (%). 


Trend 
12 months starting in January 


12 months starting in April 
3-92 3-33 
3°5 4-4 
5 5 
12 months starting in July 
05 3-35 
9 60 
2 
12 months starting in October 
3-37 


4-6 
4 


Winter half-year 
3930 337 |= 83 
3-3 28 4-0 
3 5 3 
Summer half-year 
3-74 


3-2 
4 


5 
64 
7 


Winter 


400 
1 


5 
0 


maximum entropy method to resolve into a sine wave the convex curve displayed by the (Pode Hole) 
data in Figure 1. Peaks at periodicities of 11 and 22 years, which could possibly be associated with the 
sunspot cycle, are generally absent, the closest approach being the 12 year cycle at Kew. 

5.2 Analysis of ‘mixed’ seasons. Thespectral analysis over the complete length of records was repeated 
when the data were arranged chronologically into mixed seasons as described in section 4. The larger 
number of data points than were available for the single season analysis permitted greater resolution to 
be obtained, and the power spectra were produced using 100 coefficients. As a result, more peaks were 
obtained than for the single season analysis (see section 3.3) and the larger number of data points 
available also reduced the threshold of the various levels of significance. These facts would lead one to 
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Table VIII. Significant spectral peaks of Manchester rainfall from 1786 to 1975 


Results were obtained from the Maximum Entropy Method using 60 coefficients. 
T = period (years); V = proportion of variance (%); S = level of significance (%). 

Trend ase 

12 months starting in January 
. 3-85 2:90 2-56 

3-6 ‘ 4:1 5-4 2:8 

— ° 5 4 4 
12 months starting in April 

3-87 


3-9 
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12 months starting in July 
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| 
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12 m 
4-93 
3-7 
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Winter half-year 
97 +93 
4 “7 


Summer half-year 


Winter 
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06 
4 


expect more significant peaks than were obtained with the single season analysis. The results, presented 
in Table IX show only about the same number of significant peaks as before. This indicates that fewer 
features of interest are present when the data are formed into mixed seasons than when they are arranged 
into individual seasons. The major peaks found in section 5.1 are in general absent, but another com- 
mon peak emerges at a period of around 1.14 years. This had been observed by Brunt (1927) in London 
rainfall, but only a small proportion of the variance is accounted for. 

5.3 Analysis of 80 year epochs. The single season form of spectral analysis was repeated when thecom- 
plete records were divided into 80 year epochs. This enabled one to examine the stability of the signifi- 
cant spectral peaks found in section 5.1 and, more generally, to see how the structure of the variance 





Meteorological Magazine, 108, 1979 


Table IX. Significant spectral peaks in rainfall summed over consecutive three month periods 


Results were obtained from the Maximum Entropy Method using 100 coefficients. 
T = period (years); V = proportion of variance (%); S = level of significance (%). 


Trend 
England & Wales 1727-1975 

1-15 1-00 
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evolved with time. The seasons examined were spring, summer, autumn and winter, together with the 
calendar year and the winter half-year, and the number of coefficients used was 26 (i.e. N/3). The 
reduced number of coefficients led to a smaller number of peaks being produced than for the analysis of 
the complete records, and the smaller number of data points involved also increased the thresholds of 
the various significance levels. As a result, the number of significant peaks obtained was less than in 
section 5.1. The results, presented in Tables X—XIII, show that no spectral peak achieved 5 per cent 
significance in all the epochs examined. There were also large differences in the power spectra between 
one epoch and another. This is illustrated in Figure 5 for the winter half-year rainfall over England and 
Wales. 

The differences in the structure of the variance between one epoch and another are not surprising. 
Consider a model of rainfall variance in which random fluctuations predominate, but in which regular 
oscillations are also present. Over any given frequency range, the power or variance may be regarded as 
being composed of random and regular components. The regular cyclic component may be assumed (in 
the model) to be constant with time, but the random component will vary with respect to epoch. In 19 
epochs out of 20, the random component of the variance will be contained by the lower and upper 5 per 
cent levels of significance, but over a narrow frequency band, these limits will be wide. In a series of 250 
data points, variations in the mean noise level of a power spectrum, i.e. the variance over the complete 
frequency range 0 to 0.5, will be determined by the standard error of the variance derived from 250 
points. The standard error of the variance over 1/10th of the frequency range will be equivalent to the 
standard error over the whole frequency range based on only 1/10th of the data, i.e. 25 points. The 
regular fluctuations found in section 5.1 generally occupied even narrower frequency bands. Thus the 
power over a narrow frequency band which contains a regular cyclic fluctuation is still subject to large 
random variations. When the random component of power is large, the spectral peak due to the regular 
periodicity is boosted, while when the random component is small, it is depressed. It is therefore not 
surprising that the spectral peaks found in section 5.1 did not maintain a regular level of significance 
over epochs as short as 80 years. 
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Figure 5. Maximum entropy power spectra (26 coefficients) of winter half-year rainfall over England and Wales during 
the epochs 1736-1815, 1816-95, and 1896-1975. 


Tables X—XIII show that the most regular of the spectral peaks are those at 2.1 and 2.4 years. The 
2.4 year cycle is the most important, and is most developed in the first and last epochs. The 3.9, 5 and 6 
year cycles do not receive regular support, but the 5 year periodicity in winter half-year rainfall is very 
pronounced in the most recent epoch (see Figure 4). An 11th-order filtered series of the 5-year cycle for 
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Table X. Significant spectral peaks in England and Wales rainfall during 80 year epochs 


Results were obtained from the Maximum Entropy Method using 26 coefficients. 
T = period (years); V = proportion of variance (%); S = level of significance (%). 
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England and Wales rainfall summed over 11 months starting in October is displayed in Figure 6. The 
fluctuation is seen to have been well dev.ioped between 1860 and 1885, and after 1925, since when the 
filter has passed 25 per cent of the variance. Over the complete length of the England and Wales record, 
the filter has passed 15 per cent of the variance, which is significant at the 5 per cent level. 

Another feature of the winter half-year rainfall which is unlikely to be permanent, but which has been 
well developed this century, is an 11 year periodicity. This is illustrated in Figure 7 by means of a sixth- 
order filtered series which has accounted for one-third of the total variance since 1908. The fluctuations 
have not, however, maintained a constant phase relationship with the sunspot cycle. 


6. Conclusions 


Four long-period rainfall records for England and Wales have been subjected to a maximum entropy 
spectral analysis. The data used, the spectral techniques employed and the method of calculation of 
variance and significance attributed to each fluctuation, are all imperfect and open to many criticisms. 
Unitary filters, however, were used to support the main findings of the maximum entropy analysis, 
while the effect of imperfections in the data were lessened by the fact that four independent series were 
used. 
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Table XI. Significant spectral peaks in Kew rainfall during 80 year epochs 


Results were obtained from the Maximum Entropy Method using 26 coefficients. 
T = period (years); V = proportion of variance (%); S = level of significance ( %). 
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It was found that rainfall fluctuations were largely random, and although some regular periodicities 
were evident, their use for forecasting was limited because they accounted for only a few per cent of the 
total variance. Cycles that could be regarded as permanent were located at periods of 2.1 and 2.4 years, 
of which the latter, associated with the summer half-year, was most important. Cycles that were less 
prominent, but which attained 5 per cent significance over the complete length of all the series examined, 
were found at periods of 3.9, 5, and 6 years. The 6 year cycle related mainly to the summer half-year, 
while the 5 year cycle, which has been well developed this century, is associated with the winter half- 
year. 
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Table XII. Significant spectral peaks in Pode Hole rainfall during 80 year epochs 


Results were obtained from the Maximum Entropy Method using 26 coefficients. 
T = period (years); V = proportion of variance (%); S = level of significance (%). 
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Table XIII. Significant spectral peaks in Manchester rainfall during 80 year epochs 


Results were obtained from the Maximum Entropy Method using 26 coefficients. 
T = period (years); V = proportion of variance (%); S = level of significance (%). 
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551.578.45(423+429) 
The blizzard of 18-19 February 1978 in south-west England and South Wales 


By Joyce Laing 
(Meteorological Office, Bracknell) 


Summary 


Heavy falls of snow accompanied by strong to gale force winds affected the south-west of England and South Wales 
on 15-16 February 1978 and again on 18-19 February 1978. The sudden onset and the severe drifting of the snow 
caused considerable disruption to communications and damage to property and livestock. Although some hill tops 
and open fields were scoured of snow cover, drifts of 6 metres were reported in valleys and deep-set roads. A comparison 
with past snowstorms indicates that not since 1891 has there been so much snow in conjunction with high winds over 
such a wide area of south-west Britain. 


1. Introduction 


During the third week of February 1978 snowfalls were considerable in south-west Britain and culmin- 
ated in the blizzard of 18-19 February when the combination of strong to gale force easterly winds and 
heavy snowfall practically brought life to a standstill from Hampshire to East Cornwall and north to 
Avon and Glamorgan. A quick thaw on the 20th and 21st produced some flooding on the south coast 
but in some areas drifts remained until May. 


2. Synoptic situation 


The broad pattern remained virtually unchanged over the British Isles from 15 to 21 February 1978. 
A large depression, centred at about 50°N, dominated the North Atlantic, and a ridge of high pressure 
extended south from the polar regions into central Europe, which resulted in a south-easterly airstream 
over much of the United Kingdom. 

On the 15th an occluding frontal system approached south-west England, bringing the warm, moist 
air from the Atlantic into juxtaposition with the established cold air. A small low-pressure centre 
formed on this front between the Isles of Scilly and Brest by midnight on the 15th and moved into 
northern France. A further low-pressure centre developed off south-west Ireland late on the 16th and 
also moved into northern France. Late on the 18th another occlusion moved in from the west with a 
small wave developing into a centre near Scilly by midnight. This low became stationary in the mouth 
of the English Channel on the 19th and gradually filled. Figure 1 shows the synoptic situation at 06 
GMT on 19 February 1978. 

Strong to gale force east to south-easterly winds over the south and west of the country were produced 
between the ridge which remained over northern and eastern Britain and the succession of lows near 
south-west England. On the 20th this ridge moved eastwards allowing the weakening fronts to move 
north-eastwards and the winds to moderate. 


3. Data used 


The climatological returns received in the Meteorological Office from the network of observing 
stations were the main source of information for this survey. These returns are from both official 
meteorological offices and from auxiliary and climatological co-operating stations. Use has also been 
made of the reports from the contributors to the Snow Survey of Great Britain (Meteorological 
Office, 1978). Figure 2 shows the locations of all the stations in southern England and South Wales 
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Figure 1. Synoptic situation at 06 GMT on 19 February 1978. 


which were used in this report. Supplementary accounts have also been received from the meteoro- 
logical offices at Mount Batten, Rhoose (Cardiff/Wales Airport) and Upavon, and reports in the 
national and local newspapers have been consulted. 

References to accounts of past snowstorms are given at the end of this paper. 


4. Weather 


The weather at the beginning of February 1978 was changeable with periods of mild, stormy weather 
in the south-west. From the 7th it turned gradually colder and during the three days from the 8th to the 
10th slight snow showers fell over many parts of the south-west. After a fine, sunny day on 11 February 
the next few days saw some moderate falls of snow, mainly over the higher ground. More general snow 
fell on the nights of 15/16 and 16/17 February, giving an appreciable cover over most of southern 
England, from Surrey and Sussex westward, and South Wales; more than 25cm were recorded on 
Exmoor and Dartmoor, and drifts of up to 10 ft (3 metres) quickly formed in the strong easterly winds, 
blocking many roads. 

The winds moderated on 17 February and the weather that day was fine and sunny. As a result of 
valiant efforts by council workers all main roads in the west were clear of snow by that night. However, 
after a cold night, Saturday 18 February became cloudy, and rain and sleet in Cornwall turned to snow 
as it spread further east and north. In the transitional state between rain and snow, generally in south- 
east Cornwall and south-east Devon, some unusual forms of precipitation were noted, variously 
described as ice pellets, ice needles, wet icicles or hail. The snow spread into south-west England and by 
evening had extended to South Wales and central England. The winds increased again to strong to near 
gale force from an easterly direction, causing severe drifting of the snow especially over the higher 
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ground. This blizzard with mean winds of 25-30 knots and blowing snow continued on to 19 February, 
causing havoc to communications, power supplies and livestock, mainly in Devon, Dorset, Somerset 
and Glamorgan. Temperatures were low away from the south coast; at Cardiff/Wales Airport, the 
temperature did not rise above —0.7 °C and on Exmoor the maximum recorded was —1.6 °C. South 
and west Cornwall, however, had comparatively little snow. The gale force winds continued throughout 
the 19th and the blown snow prevented accurate assessment of the amounts of snow actually falling. 
Wind speeds in gusts of 50-60 knots were recorded at a number of places and a maximum gust of 73 
knots was recorded at Gwennap Head (Cornwall). Exposed areas such as North Hessary Tor, Dart- 
moor, were scoured of lying snow while some adjacent areas had drifts of 20 feet (6 metres) or more. 
Accumulations of level snow to a depth of about 60 cm were recorded in places on Dartmoor and 
Exmoor and in Glamorgan. The greatest depth recorded was 85cm at Nettlecombe (Bird Hill) in 
north-west Somerset. 


The area affected by the blizzard stretched from the South Downs west of Beachy Head and including 
the Isle of Wight, across Dorset, Devon, Somerset and Avon and into Gwent, Glamorgan and parts of 
Dyfed. Table I gives the snow depths recorded at a selection of places throughout the area between 13 


Table I. Depths of level undrifted snow in centimetres at 09 GMT, 13-23 February 1978 


Date 13 14 15 16 17 18 19 20 21 
Altitude 
metres 
Place 
South Wales 
Aberporth 134 
Brawdy 111 
Brynamman 183 
Merthyr Tydfil 235 
Rhoose 67 


Somerset & Avon 
Filton 59 
Nettlecombe 

(Bird Hill) 96 
Hawkridge 314 
Crewkerne 101 
Yeovilton 18 


Devon 

Mount Batten 27 
Exeter 32 
Starcross 9 
Chawleigh 

Princetown 414 
Burrator 


Dorset 
Dorchester 
Winfrith 
Poole 
Hurn 


Wiltshire 
Boscombe Down 


Hampshire 
Southampton 3 
St Catherines I.0.W. 16 15 


Figures in brackets are estimated values; tr = trace, i.e. less than 0-5 cm;—indicates no observation; see Figure 2 for 
location of stations. 
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Figure 3. Maximum snow depths (level undrifted snow) measured at 09 GMT on 19 or 20 February 1978. 


and 22 February 1978 and Figure 3 shows the greatest depths recorded on the mornings of 19 and 
20 February, together with the extent of snowfall on 19 February. 

On the 20th the winds slowly moderated and further snowfall in the west turned first to freezing rain 
and then, as temperatures rose, to rain, which, together with the melting snow, brought further prob- 
lems with flooding in some areas. Although much of the snow thawed quickly, drifts remained for some 
time in places in the Glamorgan valleys and on Exmoor and Dartmoor; even as late as 1 May a drift two 
feet deep still remained on Dartmoor. 

The sequence of events as they occurred at Rhoose between 12 and 20 February is given in Figure 4. 
This shows the rise in wind speeds coinciding with the snow on 15-16 and 18-19 February and also the 
combination of below-freezing temperatures and strong winds on 18-19 February. In fact at Rhoose 
there were 15 hours on 18-19 February 1978 when the temperature was below freezing and the associ- 
ated hourly mean wind speed was 29 knots or more. This is no common occurrence and there have 
been only 30 such hours during the past 18 years; only one of these occasions was in February (February 
1969, also in conjunction with snow.) 


5. Some problems caused by the blizzard 


Eyewitnesses in some of the worst affected areas spoke of the suddenness with which vehicles became 
buried in the snow; the drifts were forming much more quickly than they could be cleared. At one stage 
all inessential transport was prohibited from entering Devon to ease the work of the snow ploughs 
and rescue workers. 

Farmers in Devon and Somerset reported drifts up to the eaves of their houses and on Dartmoor 
some drifts reached the top of the walls of Princetown Prison. Sheep, which tend to seek shelter in the 
lee of banks or hedges, were buried in the drifts, though many were rescued after several days. Un- 
fortunately the floods which followed the blizzard drowned some of the sheep before they could be dug 
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out. A farmer near Lynton had so much snow in his yard that he could not find his horses. Later he 
managed to dig down to them where they were lying on their backs and get them into the barn where 
they thawed out and revived. In a number of places milk could not be transported from the farms and 
had to be poured away. 

Considerable damage was done to trees and property. Glasshouses suffered greatly as they collapsed 
under the weight of snow. On the west side of Dartmoor hundreds of trees were brought down by the 
snow and wind. The Royal Mint, Llantrisant, had to close and the staff were sent home. Electricity 
lines were damaged and this caused pump failure and lack of water supplies to some areas. In the Cardiff 
area many people attending the Wales—Scotland rugby match found themselves stranded for several 
days. 

The thaw brought its attendant problems. Dense fog curtailed helicopter flights that were engaged in 
rescue work and in flying in supplies to stranded communities and fodder for cattlz and other livestock. 
Heavy rain combined with the melting snow coinciding with the high tides brought flooding to Kings- 
bridge (Devon) where the water rose to seven feet in the town. Fields were flooded in many areas, 
especially in the Taw and Torridge valleys and around Taunton. Huge slabs of snow drifted down the 
swollen streams, and lakes formed behind snowdrifts. 

The weather observers, both professional and amateur, had an arduous task during these days as 
instruments became buried in snow and the difficulties of reaching the site became almost insurmount- 
able. At Rhoose the enclosure was inaccessible because of deep snow from the 18th to the 20th. At 
Hartland Point on the 18th observers were unable to open the screen because of icing and remarks for 
the 19th—20th state ‘station closed down—snowed under’. The climatological observer at Sidmouth 
turned out although suffering from bronchitis but found the journey to the site too strenuous and had to 
seek shelter in a nearby house. The following account by the observer at Chawleigh (Devon) gives some 
idea of the appalling conditions. 

‘Doing the 2100 ob at Chawleigh on Saturday 18th was very difficult as it was next to impossible to see; 
however the 40 metre journey to the enclosure was made successfully. The 0900 [observation] on Sun- 
day morning however was done in impossible conditions—great drifts of snow 15 ft high or more had to 
be climbed with visibility almost non-existent, the temperature was —2.0 °C and the wind near 50 kt— 
indeed conditions were almost beyond the limit of human endurance and the observation was a night- 
mare, as was the 2100 the same day.’ 


6. Past snowstorms in the south-west 


The following are some of the more severe snowstorms to which the south-west of England and South 
Wales have been subjected. 
1962/63. Although the total amount of precipitation during this winter was not great there were some 
heavy falls of snow in the south-west. Tredegar (Gwent) recorded 65 inches (156 cm) of level snow on 8 
and 9 February. This weather was caused by troughs of low pressure moving into southern England 
from France (Shellard, 1968; Meteorological Office, 1963). 
1947. A trough moving west from the continent covered virtually all Britain with snow on 27/28 
January, and a depression which formed on this trough became centred in the western Channel and 
prolonged the snow in south-west England, including the Isles of Scilly, where 7 inches (18 cm) were 
measured (Douglas, 1947; Meteorological Office, 1947). 
1945. A complex situation of small secondary lows off the south-west coasts moved east and brought 
snow to the Bristol Channel areas during the period 22-25 January. Cardiff recorded 30 inches (76 cm) 
of level snow; Okehampton (Devon) reported a greater depth than at any time since 1895 but remarked 
that there was little drifting (Bonacina, 1950; Jackson, 1977). 
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1933. Snowfall, associated with a polar low which developed in St George’s Channel and moved south- 
wards, affected a wide area of England and Wales on 23-26 February. In South Wales the general depth 
of level snow was 1-2 feet (30-60 cm) with deep drifts (Bonacina, 1937). 
1929. On 16 February heavy snowstorms occurred in west England and Wales. Cardiff had great 
difficulty in keeping the streets clear of snow. On Dartmoor it was estimated that up to six feet (180 cm) 
of snow fell in 15 hours. This is thought to be the deepest fall of snow in a single storm anywhere in the 
British Isles at so low an elevation as 1000 ft (300 m), but it cannot be classed as a blizzard as there was 
little wind at the time (Bonacina, 1937). 
1927. Heavy snow over the Christmas period was accompanied by gale force north-easterly winds 
engendered by a depression in the western Channel. Drifts of 20 ft (6 m) were reported on Salisbury 
Plain and on Dartmoor the prison was unapproachable for days. On the south coast the effects of 
drifting were most striking; there were reports of snow being blown off the Isle of Wight over the sea. 
In terms of drifting this storm was comparable with those of 1881 and 1891 (Bonacina, 1937; Jackson, 
1977; Douglas, 1928). 
1891. From 9 to 13 March the ‘great west country blizzard’ affected the whole of the south of England 
and South Wales with areas of great intensity in Kent and in Devon and Cornwall. The average depth 
of snow in the south-west was reported as about 2 ft (60cm) with immense drifts, and locally on Dart- 
moor this amount was doubled by a secondary outbreak of snow on 12 March. Great disruption to 
transport was caused, with many trains marooned and many ships sunk in the Channel. In those days 
drifts had to be cleared by manual labour, the railway gangers vying with each other to keep their own 
lines open. This blizzard was the result of waves developing on a cold front: one such low-pressure 
centre, on 9-10 March, deepened at Ushant and moved up the Channel, depositing snow in drifts with 
the strong east to north-east winds (Bonacina, 1928; Jackson, 1977; Symons, 1891; Shellard, 1968). 
1881. The great blizzard of 18-21 January mainly affected the Wessex area, that is to say east Devon, 
Somerset, Dorset, Wiltshire and Hampshire. The Isle of Wight was particularly badly hit: ‘the depth 
of snow which was doubled on 20th was about 3 feet with gigantic drifts’. A depression was centred 
near the Channel Islands on the 18th and moved steadily up Channel. The blizzard was described at the 
time as the worst since 1836 (Bonacina, 1928; Symons, 1891). 

Figures 5 and 6 show the distribution of the maximum snowfall in 1881 and 1891 (Jackson, 1976). 


7. Conclusion 


The main features of the blizzard of February 1978 were the comparative shortness of the cold spell, 
the persistence of the snowfall in one area for two days, the severity of the drifting and the quickness of 
onset of both the snow and the thaw. 

The snowfall was the heaviest to occur in the south-west for many years. The combination of so much 
snow with gales had not been experienced in Devon since 1927 and in south-west England and South 
Wales as a whole since 1891. 

Of the eight outstanding storms in the south-west since 1850 the visitations of March 1891 and January 
1881 seem to have been very similar in extent and severity to those of 1978. All three were the result of 
low-pressure centres in the western Channel. 

Bonacina (1928) remarks that in the storms of 1881, 1891 and 1927 there were two phases to the snow- 
fall with a gap of one or two days between, and although in each phase the depths were substantial it 
was the combination of the two falls which produced such outstanding conditions. This feature was 
present also in 1978; the snow on 15 and 16 February was followed by the main blizzard of 18-19 
February which more than doubled the depths of snow and drifts in some areas. 
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Figure 5. Maximum snow depths 17-21 January 1881. 
(N.B. Correction: the ‘15’ and ‘40’ labels should be interchanged.) 
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Figure 6. Maximum snow depths 9-13 March 1891. 
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Notes and news 
Dr P. M. A. Bourke 


Dr P. M. Austin Bourke, Director of the Irish Meteorological Service since 1964, retired in May of 
last year. Apart from a few years as a university lecturer in mathematics, his career was spent entirely in 
the Irish Meteorological Service. He became well known for his work on agricultural meteorology and 
had regular close contact with Mr L. P. Smith and the Agricultural Meteorology section of the Meteoro- 
logical Office; he preceded Mr Smith as president of the WMO Commission for Agricultural Meteoro- 
logy, serving from 1958 to 1962. 

Dr Bourke is highly regarded, not only for his intellectual ability, but also for his humour and charm 
which have enlivened many international meetings. We wish him well in his doubtless active retire- 
ment. 


Obituary 


We record with regret the death on 12 January 1979 of Miss C. A. Parkhouse, Assistant Scientific 
Officer, HQ Strike Command, after a brief illness. Miss Parkhouse joined the Office in August 1978. 


Corrections 


Meteorological Magazine, Volume 107, December 1978; article by J. P. Cowley. The unit used for 
Figure 5, pp. 366-372, is MJ m-?. 

Meteorological Magazine, Volume 108, January 1979; article by C. J. Richards. In the penultimate 
line of the caption to Figure 2, p. 14, 1970 should read 1976. In equation (10), p. 22, a fraction bar and a 
figure 2 have been omitted after the second square-root sign. 
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